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Abstract 

We present the analysis of exclusive - ^0~ semileptonic heavy meson decays 
using the constituent quark model based on the light-front quantization. Our 
model is constrained by the variational principle for the well-known linear plus 
Coulomb interaction motivated by QCD. Our method of analytic continuation to 
obtain the weak form factors avoids the difficulty associated with the contribution 
from the nonvalence quark-antiquark pair creation. Our numerical results for the 
heavy-to-heavy and heavy-to-light meson decays are in a good agreement with 
the available experimental data and the lattice QCD results. In addition, our 
model predicts the two unmeasured mass spectra of 1 5'o(65) and 3 Si(bs) systems 
as M fe5 =9657 MeV and M bs = 5424 MeV. 
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In recent years, the exclusive semileptonic decay processes generated a great ex- 
citement not only in measuring the most accurate values of the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix elements but also in testing diverse theoretical approaches to 
describe the internal structure of hadrons. Especially, due to the anticipated abun- 
dance of accurate experimental data from the 5-factories (e.g. HERA-B at HERA, 
BaBar at SLAC and Belle at KEK), the heavy-to-heavy and heavy-to-light meson 
decays such as 5— >7r, D—>tt(K) etc. become invaluable processes deserving 
thorough analysis. While the available experimental data of heavy meson branching 
ratios have still rather large uncertainties^], various theoretical methods have been 
applied to calculate the weak decay processes, e.g., lattice QCD0, |3], ||, QCD sum 
rules||, Heavy quark effective theory||, and quark models]?], |8|, |, [TO], [LI], [12], [13| . In 
particular, the weak transition form factors determined by the lattice QCD[|J provided 
a useful guidance for the model building of hadrons, making definitive tests on existing 
models, even though the current error bars in the lattice data are yet too large to pin 
down the best phenomenological model of hadrons. These weak form factors, however, 
are the essential informations of the strongly interacting quark/gluon structure inside 
hadrons and thus it is very important to analyze these processes with the viable model 
that has been very successful in analyzing other processes. 

In this letter, we report the analysis of exclusive semileptonic decays of a heavy 
pseudoscalar meson into another heavy or light pseudoscalar meson using the light-front 
quark model (LFQM) which has been quite successful in the analysis of electromagnetic 
form factors, radiative decays and K—>7i transition form factors |7j, ||. The LFQM takes 
advantage of the equal LF time {r=t + z/c) quantization [[H] and includes important 
relativistic effects in the hadronic wave functions. The distinguished feature of the 
LF equal-r quantization compared to the ordinary equal-t quantization is the rational 
energy- momentum dispersion relation [HJ which leads to the suppression of vacuum 
fluctuations with the decoupling of complicated zero modes |TB], [T7J and the conversion 
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of the dynamical problem from boost to rotation |18]. The recent lattice QCD results |T9 



indicated that the mass difference between rf and pseudoscalar octet mesons due to 
the complicated nontrivial vacuum effect increases (or decreases) as the quark mass m q 
decreases (or increases), i.e., the effect of the topological charge contribution should 
be small as m q increases. This supported in building the constituent quark model @ 
in the LF quantization approach because the complicated nontrivial vacuum effect 
in QCD can be traded off by the rather large contstituent quark masses. We have 
also circumvented the problem of assigning the dynamical quantum numbers J PC to 
hadrons by using the Melosh transformation of each constituent from equal t to equal 

Moreover, one of the most distinctive advantages in the LFQM has been the utility 
of the well-established Drell-Yan-West (q + =q° + g 3 =0) frame for the calculation of 
various form factors pD| . By taking the "good" components of the current (j + and j±), 
one can get rid of the zero mode|16[ problem and compute the full theoretical prediction 
for the spacelike form factors in q + =0 frame. The weak transition form factors that we 
are considering, however, are the timelike q 2 > observables. Our method is to rely on 
the analytic continuation from the spacelike region to the timelike region calculating the 
"good" components of the current in the q + =0 frame||. If we were to take the q + ^0 
frame, then we must take into account the higher Fock-state (nonvalence) contributions 
arising from quark- ant iquark pair creation (so called "Z-graph" ) as well as the valence 
configuations. In fact, we notice that a few previous analyses[Il]] were performed in 
the q + y^0 frame without taking into account the nonvalence contributions. We find 
that such omission leads to a large deviation from the full results ||. Our method is 
to rely on the analytic continuation from the spacelike region to the timelike region 
calculating the "good" components of the current in the q + =0 frame. 

The key idea in our LFQM[]7| for mesons is to treat the radial wave function as a trial 
function for the variational principle to the QCD-motivated Hamiltonian saturating 
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the Fock state expansion by the constituent quark and antiquark. The spin-orbit wave 
function is uniquely determined by the Melosh transformation. We take the QCD- 
motivated effective Hamiltonian as the well-known linear plus Coulomb interaction 
given by 

H qq = H + V qq = yfm* + k 2 + ifiq + k 2 + V qq , (1) 

where 

V q ^V + V^a + W-f r + %^^. (2, 

We take the Gaussian radial wave function <j)(k 2 ) = iVexp(— k 2 /2/3 2 ) as our trial wave 
function^ to minimize the central Hamiltonian J7|. Among the light-quark mass and 
the potential parameters {m u {= md),/3 u d(= Puu),o,,b, k}, only 4 parameters are in- 
dependent because of the constraint from the variational principle. Furthermore, the 
string tension 6=0.18 GeV 2 and the constituent u and d quark masses m u =md=0.22 
GeV are rather well known from other quark model analyses commensurate with Regge 
phenomenology|S|]. Thus, using the experimental values of p and n masses and the vari- 
ational constraint, we can fix the remaining parameters a,n, and /3 n j as a = —0.724 
GeV, (3 ud - = 0.3695 GeV, and k = 0.313, respectively JTp. More detailed procedure of 
determining the model parameters of light-quark sector (u and s) can be found in 0. 
It is very important to note that all other model parameters such as m c , m^,, (3 UC , j3 u b, 
etc. are then uniquely determined by the same procedure as the light-quark analysis 
0. The procedure of determining model parameters constrained by the variational 
principle is shown in Fig. 1, where the lines of qq and qc (q=u and d) etc. represent 
the sets of {m q ,m q , (3 qq \ and {m q , m c , /3 qc }, respectively, etc. Because all the lines in 

1 Even though one can in principle expand the radial function 4> n .i=o(k 2 ) with a truncated set of 
HO basis states our choice of radial wave function turns out to be sufficient for the analysis of the 
ground state ^ and 1 ground state meson properties |t]]. 
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Fig. 1 should go through the same point of (6=0.18 GeV 2 ,/t = 0.313), the parameters 
of m c , nib, /3 UC , (3 u d, etc. are all automatically determined without any adjustment. 
Our model parameters obtained by the variational principle are summarized in Table 
1. 

Our predictions of the ground state meson mass spectra and the decay constants of 
various heavy pseudoscalar mesons are summarized in Tables 2 and 3, respectively, and 
compared with the available experimental datafll]] and the lattice QCD result s@]. Our 
predictions of ground state meson mass spectra agree with the experimental data|3J 
within 6% error. Furthermore, our model predicts the two unmeasured mass spectra of 
1 S {bb) and 3 Sx(bs) systems as M fe5 =9657 MeV and M b - s =hA2A MeV, respectively. Our 
values of the decay constants are also in a good agreement with the results of lattice 
QCD[[| anticipating future accurate experimental data. 

The matrix element of the current j M = q 2 ^^Q\ for 0~(Qig)^0~(g29) decay is given 
by two weak form factors / + and /_, viz., 



where = (Pi — P^)^ is the four-momentum transfer to the lepton and m 2 < q 2 < 
(Mi — M 2 ) 2 . In the heavy quark limit Mi( 2 )— >oo, the form factor f+(q 2 ) is reduced 
to the universal Isgur-Wise (IW) function, £(ui • v 2 ) = [2^M 1 M 2 /(M 1 + M 2 )]f + (q 2 ), 
where V\^) — Pi(2)/Mi^)- In LFQM, the matrix element of the weak vector current 
can be obtained by the convolution of initial and final LF meson wave functions in 
q + =0 frame: 



(P^QilPi) = U{q 2 ){Pi + P2Y + f-{q 2 )q 



(3) 



(P 2 |g 2 7 M Qi|A) 




rrii — rriq) 2 



xTr 75(^2 + m 2 )Y(A +^1)75(^5-^9) , 



(4) 
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where M 2 = (k\ + mf)/(l — x) + (£r[ + In Eq. (4), the form factors f+{q 2 ) 

and f-(q 2 ) are obtained by taking the "good" components of the current (j + and j±). 
The detailed derivation of f±(q 2 ) can be found in |J. Our analytic continuation to the 



timelike region has verified H the equivalence to the dispersion method [[12 



Our numerical results of the decay rates for D^tt(K), D s ^f](i] f ), and B — > tt(D) 
processes are consistent with the experimental data as summarized in Table 4. It is 
interesting to note that our value of rj-r]' mixing angle, 0su(3)=— 19° presented in 0, are 
also in agreement with the data for D s ^r](r]') decays. One should note that the number 
of events for the D — > ir data is currently very small compared to other pro cesses 0. 

In Figs. 2(a) and 2(b), we present the form factors f+ K {q 2 ) and fo' K {q 2 ), re- 
spectively, with the definition of f (q 2 ) as f (q 2 ) = f + (q 2 ) + q 2 f-(q 2 )/(M? - Mf), 
and compare with the experimental data as well as the lattice QCD results^]. Note 
that our value of f+ K (0) =0.736 is within the error bar of the measured value 0, 
_/+ xpt '(0) = 0.7 ± 0.1. In Fig. 3, we present the form factor f+ n (q 2 ) and compare 
with the results from lattice QCD [[]]]. Our result is very close to the UKQCD[[| re- 
sults for a wide range of momentum transfer. In Fig. 4, our prediction of the IW 
function for B^-D transition are compared with the experimental data[]2"T|, P^ ]. Our 
prediction of the slope p 2 =0.8 of the IW function at the zero-recoil point defined as 
£(vi ■ V2) = 1 — p 2 {vi ■ v 2 — 1) is quite comparable with the current world average 
Pavg.=0.66±0.19[jll extracted from exclusive semileptonic B^>Dtv decay. 

In conclusion, in this paper, we analyzed the exclusive 0~— >0~ semileptonic heavy 
meson decays using the LFQM constrained by the variational principle for the QCD- 
motivated effective Hamiltonian with the well-known linear plus Coulomb interaction. 
Our model not only provided overall a good agreement with the available experimental 
data and the lattice QCD results for the weak transition form factors and branching 
ratios of the heavy-to-heavy and heavy-to-light meson decays but also rendered a large 
number of predictions to the heavy meson mass spectra and decay constants. Our 
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model in fact predicted the masses of heavy mesons, i.e., M b 5( 1 S'o)=9657 MeV and 
Mbs( 3 Si)= 5424 MeV. We have overcome the difficulty associated with the nonvalence 
Z-graph contribution in timelike region by the analytic continuation of weak form fac- 
tors from the spacelike region. Our numerical computation confirmed the equivalence 
of our analytic continuation method and the dispersion relation method Hl2|pi We think 



that the success of our model hinges on the advantage of light-front quantization re- 
alized by the rational energy-momentum dispersion relation. It is crucial to calculate 
the "good" components of the current in the reference frame which deletes the com- 
plication from the nonvalence Z-graph contribution. We anticipate further stringent 
tests of our model with more accurate data from future experiments and lattice QCD 
calculations. 
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Research Scientific Computer Center are also acknowledged for the grant of supercom- 
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Table 1: The constituent quark masses m[GeV] and the Gaussian parameters /3[GeV] 
for the the linear potential obtained by the variational principle. q=u and d. 



m q 


m s 


m c 


m b 




As 


PqS 


Pq-C 


0sc 


Pec 


Pqb 


Psb 


Pbb 


0.22 


0.45 


1.8 


5.2 


0.3659 


0.4128 


0.3886 


0.4679 


0.5016 


0.6509 


0.5266 


0.5712 


1.1452 



Table 2: Fit of the ground state meson masses [MeV] with the parameters given in 
Table I. Underline masses are input data. The masses of (o> — 0) and (tj — t/) were used 
to determine the mixing angles of u — <fi and rj — rj' [7] , respectively. 



'So 


Expt.i 


Prediction 


3 ^ 


Expt.i 


Prediction 


71 


135 ±0.0006 


135 


P 


770 ± 0.8 


770 


K 


498 ± 0.016 


478 


K* 


892 ± 0.26 


850 


V 


547 ± 0.12 


547 




782 ± 0.12 


782 


v( 


958 ±0.14 


958 





1020±0.008 


1020 


D 


1865±0.5 


1836 


D* 


2007± 0.5 


1998 


D s 


1969±0.6 


2011 


D* 


2112±0.7 


2109 


Vc 


2980±2.1 


3171 


J/4 


3097±0.04 


3225 


B 


5279± 1.8 


5235 


B* 


5325± 1.8 


5315 


B s 


5369±2.0 


5375 


(bs) 




5424 


(65) 




9657 


T 


9460± 0.21 


9691 
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Table 3: Decay constants [Me V] for various heavy pseudoscalar mesons. 



References 


fo 


fD s 


fa 


/b s 


Ours 


139.2 


164.8 


121.2 


144.2 


Lattice |2 


1 


141.4±21.2 


155.6±21.2 


120.2±24.7 


137.9±24.7 


Expt.(T[ 




< 219 


137 - 304 







Table 4: Form factors /+(0) and branching ratios (Br.) for various heavy me- 
son semileptonic decays for 0~^0~. We use 0g^L=— 19° for D s — > r)(r]') de- 
cays and the following CKM matrix element: | V cs |=1.04±0.16, | V cd |=0.224±0.016, 
|K&| = (3.3±0.4±0.7)xlQ- 3 , and | Vj )C |=0.0395±0.003 [1], 



Processes 


/+(o) 


Br. 


Expt.i 


D^K 


0.736 


(3.75 ± 1.16)% 


(3.66 ±0.18)% 


D -> 7T 


0.618 


(2.36 ±0.34) x 10~ 3 


(3.9i^±0.4) x 10~ 3 




0.421 


(1.8 ±0.6)% 


(2.5 ±0.7)% 




0.585 


(9.3 ±2.9) x 10~ 3 


(8.8 ±3.4) x 10" 3 


B ^ 7T 


0.273 


(1.40 ±0.34) x 10~ 4 


(1.8 ±0.6) x 10" 4 


B -> D 


0.709 


(2.28 ±0.20)% 


(2.00 ±0.25)% 
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Figure 1: The parameters m s , m c , rrib, (3 qs , (3 qc , etc. satisfying variational principle. The 
qq and qc etc. represents the sets of (m g ,m g , j3 qq ) and {m qi m C) j3 qc ) etc., respectively. 
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Figure 2.b: The form factor fo(q 2 ) for D^K transition compared with the lattice QCD 
results [4]. 
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Figure 3: The form factor f+{q 2 ) for B^tt transition compared with the lattice QCD 
results [3]. 



15 



0.15 




Figure 4: The IW function £(vi ■ v 2 ) for B^D transition compared with the experi- 
mental data of ARGUS [21] (square) and CLEO [22] (circle). 
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